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ABSTRACT 


The use of constrained viscoelastic layer damping treatments has 
recently been found to be effective in reducing unwanted resonant 
vibrations in structures. The proper design of such damping treatments 
has been difficult, however, since closed form equations to describe system 
behavior can only be used for simple beams and plates. In addition, finite 
element analysis has been difficult since the material properties of 
viscoelastic materials show a marked dependence on frequency and 
temperature, resulting in the need to run an excessive number of finite 
element analyses for each design. Two recently developed design tools are 
used in this paper to design a proper constrained viscoelastic layer 
damping treatment for a test structure consisting of two concentric 
cylinders connected by four vanes. The first design tool is subsystem 
analysis, where a portion of the test structure is approximated by 
appropriate closed form equations to determine system behavior with 
different applied damping treatments. The second design tool is a finite 
element technique called the Modal Strain Energy method. 

After application of the constrained viscoelastic layer damping 
treatment to the outer shell of the test structure, vibration response 
amplitudes were found to drop significantly for all frequencies except in the 
lowest fifth of the frequency range. This was due to the low shear modulus 
of the viscoelastic material at these low frequencies, which made the 
viscoelastic material unable to provide good shear resistance. Structural 
damping also increased. Further application of damping treatments to the 
vanes of the structure produced only slightly more damping and response 
amplitude reduction. 

The Modal Strain Energy method did not show good agreement with 
the experimental results. In general, the Modal Strain Energy method was 
a conservative design tool in that it predicted higher responses and less 
damping for the test structure than was obtained experimentally. Errors in 


1V 


the Modal Strain Energy analysis were due to coarseness of the finite 
element model and inaccuracy in the modeling of the bolted joints of the 
structure. 
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I. INTRODUCTION 


A. BACKGROUND 

Modern high strength materials such as steel, aluminum and 
titanium are widely used in all areas of construction. However, these 
materials have very little inherent damping. Structures constructed of 
these materials are apt to experience large, unwanted resonant vibrations 
which can result in structural damage and the transmission of large 
amounts of noise to the environment. Most materials which have high 
inherent damping, such as viscoelastic materials, have very low load 
carrying capability, making them unsuitable for present day high strength 
applications. A combination of a high strength material with a material 
that has high inherent damping can exhibit the necessary 
strength/damping characteristics, if properly designed. One of the more 
effective methods of combining these materials is through the use of 
constrained viscoelastic layer damping treatments. Unfortunately, the 
proper design of constrained viscoelastic layer damping treatments for 
structures can be very difficult for a number of reasons, including the 
following: 


e The material properties of viscoelastic materials vary greatly with 
temperature and frequency. 

e Closed form equations which allow easy analytical solutions to 
design problems are only available for simple beams and plates. 

e The amount of damping that can be obtained from a given 
constrained viscoelastic layer damping treatment is heavily 
dependent on the location of the treatment on the structure. 
Improper placement of the treatment on the structure will greatly 
mitigate the potential improvements in system damping. 

e Until recently, the variation of viscoelastic material properties with 
frequency made finite element modeling a prohibitively costly way of 


determining structural system damping due to the large CPU time 
requirements. This was especially so if an iterative design process 
was to be used, requiring analyses of many damping treatment 
combinations. 


BB: PURPOSE 

Recently, a new finite element analysis technique has been developed 
which is effective in modeling structures which use constrained 
viscoelastic layer damping treatments, and is cost effective enough to use 
with complicated structures. This technique, called the Modal Strain 
Energy (MSE) method, uses structural strain energy to approximate the 
damping properties of the structure with an applied damping treatment. 
The MSE method has been validated by comparison to classical frequency 
response finite element techniques and by experimentation for many 
simple structural geometries [Refs. 1,2]. The purposes of this study are to 
develop and apply a design procedure for constrained layer damping 
treatments for general, complicated plate and shell type structures, and to 
investigate the effectiveness of a constrained layer damping treatment for a 
typical plate and shell type structure. Experimental techniques are used to 
measure the vibrational amplitudes in the structure to determine the 
effectiveness of the constrained layer damping treatment, and to compare to 
the results obtained using the MSE method. 


IL, THEORY 


A. EQUATIONS OF MOTION AND VIBRATION MODES 
Continuous systems that exhibit only one principal vibrational mode 
(such as a spring and viscous damper with an attached mass) are said to 


have only one degree of freedom. The equation of motion for such systems is 
as follows: [Ref. 3] 


mx(t)+ cx(t)+ kx(t) = F(t) (Zap 

where, 

x(t) = system displacement 

x(t) = system velocity 

x(t) = system acceleration 

m =system mass 

c =system damping 

k =system stiffness 

F(t) = externally applied time varying force 


The external force and resulting response can be assumed to be of the form: 


F(t) = Fe! (2.2) 
x(t) = XeI* (2.3) 
where, 
F = force amplitude 
X = response amplitude 
@ = frequency of harmonic load (radians/sec) 


Substitution of the assumed force and response into the system equations of 
motion provides the solution for the displacement X: 


= = a 2) , (2.4) 
—-mo’) + j(ca) 

If the continuous system cannot be represented by only one degree of 
freedom, two options are available to describe the system's vibratory 
behavior. Closed form solutions to the equations of motion can be found for 
simple geometries such as beams and plates. If the system is more 
complicated than that, however, the system response can only be found by 
discretization of the equations of motion. In this case, the system is 
physically split into many small pieces whose motion is then solved for. The 
equations of motion for this case are as follows: 


[M] {x} +[(C] {x} + [K] {x} = {F(t)} (2.5) 

where, 

{x} = system displacement vector 

{x} = system velocity vector 

{x} = system acceleration vector 

[M] = system mass matrix 

[C] = system damping matrix 

[K] = system stiffness matrix 

{F(t)} =externally applied load vector 


The solution to this discretized form of the system equations of motion can 
be calculated by hand if the system can be adequately described by two or 
three degrees of freedom. For more complicated systems, however, the 
matrix operations become too complicated and must be solved using a 
digital computer. 

Modal analysis involves the determination of system natural 
frequencies and mode shapes of vibration for each natural frequency. The 
natural frequencies and mode shapes are solved for as follows. First, the 


applied force vector and damping matrix are set equal to zero. This leaves 
the equations of motion in the form: 


[M] {x} +[K] {x} = 0 (2.6) 


This is a simple eigenvalue problem that is easily solved for the eigenvalues 
and eigenvectors. For the physical system being modeled, the eigenvalues 
will represent the system's natural frequencies, while the eigenvectors will 
represent the system's mode shapes. The system's equations of motion can 
be transformed into a set of uncoupled equations which can then be solved 
individually through the use of the following transformation equation: 


ix(t)} =[o] 1a(t)} (2.7) 


where, 
{a(t)}= principal or modal coordinate vector 
[>] = mode shape matrix 
{x(t)} = original displacement vector 


Assuming that the damping matrix in equation (2.5) can be approximated 
through the use of an appropriate diagonal matrix, and using the following 
two relations: 


or =< (2.8) 
( = — 
2+/km (229) 


a new system of individual uncoupled equations can be obtained [Ref. 2]: 


&(t) + 20.0, 6.(t) + o%a(t) = F(t) (2.10) 
where, 


Oo .(t) = modal acceleration for uu mode 
a .(t) = modal velocity for ead 

a.(t) = modal displacement for erode 
¢, =critical damping ratio for ag mode 
@.= igs natural frequency 


= h 
f .(t) = modal force for A mode 


Since this system of equations is uncoupled, the equations are readily solved 


one at a time. 


B. VISCOELASTIC DAMPING 

One of the effective means of controlling unwanted resonant 
vibrations is through the use of viscoelastic materials. Many different 
applications of viscoelastic damping have been developed including free 
layer viscoelastic damping treatments, constrained layer viscoelastic 
damping treatments, concentrated dampers, and others. Every specific 
application takes advantage of the viscoelastic material's unique ability to 
dissipate vibrational energy as heat through the stretching and relaxation 
of long polymeric molecule chains [Ref. 4]. Since the material dissipates 
energy through the stretching and relaxation of its molecules, anything 
that affects the ability of the molecules to move will affect the damping 
capability of the material. Thus, the material properties of viscoelastic 
materials, especially shear modulus and loss factor, vary considerably with 
vibration frequency and temperature. The shear modulus is one third of the 
storage (Young's) modulus, while loss factor is defined as twice the critical 


damping ratio. A plot of storage modulus and loss factor vs frequency for 


viscoelastic materials is shown in Figure 2.1 [Ref. 4]. 


Real part of the modutus, E 
Loss factor n 


Frequency (log scale) 


Figure 2.1. Variation of viscoelastic material 
properties with frequency [Ref. 4]. 
Note that the storage modulus rises with increasing frequency, while 
loss factor rises and then falls. A plot of storage modulus and loss factor vs 
temperature for viscoelastic materials is shown in Figure 2.2 [Ref. 4]. 
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Figure 2.2. Variation of viscoelastic material 
properties with temperature [Ref. 4]. 


Note that storage modulus generally drops as temperature increases, 
while the loss factor rises and falls. The transition region is the most 
important region for most damping applications, since loss factor peaks 
and the storage modulus is high enough to allow good resistance to 
shearing. 

One of the more effective methods of using viscoelastic materials is 
through the use of a constrained viscoelastic layer damping treatment. In 
this application, the system to be damped is covered, either totally or 
partially, with a layer of viscoelastic material followed by a layer of some 
constraining material, as shown in Figure 2.3. 


CONSTRAINING LAYER 
VISCOELASTIC LAYER 





BASE STRUCTURE 


Figure 2.3. Constrained layer damping treatment. 


As the structure deforms the constraining layer prevents movement 
of the upper portion of the viscoelastic material. At the same time, the lower 
portion of the viscoelastic material is forced to move with the base structure. 
Thus, the viscoelastic layer is placed into shear, with the shearing action 
stretching the molecules of the viscoelastic material and causing 
dissipation of vibrational energy. Since the shearing action in the 
viscoelastic material is so great, bending deformations in the viscoelastic 


material layer can be neglected. In this case, the energy dissipation in the 
viscoelastic layer is given by the following expression [Ref. 4]: 


Diet | nNEydv (2.11) 
ol 


Most significantly, the amount of energy dissipated in the viscoelastic 
material layer is proportional to the square of the shear strain. Thus, the 
location of the damping treatment on the structure is the single most 
important factor in designing an appropriate constrained viscoelastic layer 
damping treatment. Also, of course, the energy dissipation is dependent on 
the material properties of the viscoelastic material, which themselves are 
dependent on frequency and temperature. Finally, since the energy 
dissipation integral is taken over the viscoelastic material volume, the area 
and thickness of the constrained layer damping treatment will influence 
the amount of energy dissipation realized. 


C. MODAL STRAIN ENERGY METHOD 

Most complicated structures can be adequately analyzed through the 
use of finite element techniques. Unfortunately, finding the frequency 
response of a structure whose material properties change with frequency 
has been extremely time consuming and costly, since many individual 
runs must be made with the material properties changing at each 
frequency increment [Ref. 1]. Add to this the necessity to try different design 
options depending on the specific design requirements, and finite element 
analysis can be too great of an expense. The recent development of the 
Modal Strain Energy (MSE) method, however, makes finite element 
analysis for even the most complicated structures an analysis and design 
possibility [Ref. 5). 

Central to the MSE method is the idea that the structure's damped 
mode shapes and natural frequencies can be adequately approximated by 
the undamped mode shapes and natural frequencies of the structure [Ref. 
1]. Thus, costly complex eigenvalue analysis is completely avoided. The first 
step is to find the mode shapes and natural frequencies of the structure by 


conducting a finite element modal analysis. Also, the amount of strain 
energy in each element of the finite element model is determined. It can be 
shown [Ref. 1] that the modal loss factor can be approximated by the 
following expression: 





7 Maes (2.12) 
Nig av 


where, 
i = mode number 
1; = modal loss factor 
TN, ; = viscoelastic material loss factor 
? 


V,, ; = strain energy in viscoelastic material for ith mode 
? 


h mode 


Vm ; = strain energy in total structure for it 
Since n=2¢ this expression can be readily substituted into equation (2.10). 
Also, since the modal analysis gives the system's natural frequencies and 
mode shapes, all of the constant values of equation (2.10) have been found, 
and the equations of motion can now be solved. 
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Ill. DAMPING DESIGN 


A. TEST STRUCTURE 

The test structure used for this design procedure was a set of two 
concentric cylindrical shells connected by four vanes as shown in Figure 
3.1. This structure had previously been used to investigate the ability of 
bolted connections to dissipate vibrational energy and thereby increase 
structural damping [Refs. 6,7]. This structure was useful because it was a 
structure typical of those with very little inherent damping. A constrained 
layer damping treatment was chosen for application in an attempt to 
reduce the large resonant vibrations that this structure experienced [Ref. 
6,7]. The following variables were considered in designing the proper 
constrained layer damping treatment: 


¢ Damping treatment location on the structure. 

e Temperature and frequency range of interest. 

¢ Material type for the viscoelastic and constraining layers. 

e Area and thickness of the viscoelastic and constraining layers. 


B. DAMPING TREATMENT LOCATION 

The first step in designing a proper damping treatment was to 
determine where to place it on the structure. From equation (2.11), it was 
clear that locations that would provide high shear strains in the damping 
material would be the best candidates for a damping treatment application. 

The first method of determining which areas would provide high 
shear strains in the damping material involved visual observation of the 
mode shapes of the test structure. Iverson experimentally determined the 
first eight mode shapes of vibration during his joint damping investigation 
using a Hewlett Packard 5451C Fourier Analyzer [Ref. 6]. Figures 3.2 and 
3.3 show the first two mode shapes of the eight that were determined. Table 
3.1 shows the results of this visual mode shape investigation, with a one 
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Figure 3.1. The test structure. 
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Figure 3.2. First mode shape of vibration of the test structure. [Ref. 6] 





Figure 3.3. Second mode shape of vibration of the test structure. [Ref. 6] 
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denoting a large amount of motion in a given area, and a zero denoting 
little or no motion in that area for each mode. As Table 3.1 shows, 


TABLE 3.1. VISUAL MODE SHAPE INVESTIGATION. (A "1" 
indicates a large amount of motion for that location for 
that mode. A "0" indicates little or no motion.) 





most of the structural motion took place in the outer shell area, with a 
moderate amount of motion taking place in the vanes. Almost no motion 
occurred at the vane connections to the outer or inner shells. Visual! 
observation of the animated mode shapes as seen on the HP 1347A display 
terminal reinforced these conclusions. Almost every mode analyzed 
showed sizeable deformations in the outer shell. Some modes showed 
almost no vane deformation, such as modes 1,4,7, and 8. In summary, 
investigation of the mode shapes indicated that the outer shell was a 
primary candidate for damping treatment application, with the vanes as a 
possible secondary candidate. 

The second method of determining damping treatment location was 
through the use of strain energy distribution data [Ref. 1]. In its simplest 
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definition, strain energy is the work required to deform an elastic body. For 
a rod in tension this is [Ref. 8]: 


i= SVEe, (3.1) 


XX 
where, 
U = strain energy 


EK = Young's Modulus 


Eyy = elastic strain 


V = volume 


Thus, strain energy, like the energy dissipation of equation (2.11), is 
proportional to the square of the deformation of the object in question. Thus, 
any area of the test structure that had a high strain energy was a likely 
candidate for a constrained layer damping treatment since bending 
deformations in the structure would lead to shear deformations in the 
viscoelastic material due to the constraining layer. To determine the strain 
energy distribution in the structure, a finite element model of the 
undamped structure was constructed using MSC/NASTRAN. To model 
this structure adequately, a finite element model consisting of 128 grid 
points and 108 shell elements was constructed as shown in Figures 3.4 and 
3.5. CQUAD4 shell elements were used to model the structure. A 
cylindrical coordinate system was used instead of the default rectangular 
coordinate system to make grid point location definition easier. The grid 
points were numbered as shown. A modal analysis was run, with strain 
energy for each element for each mode selected as output. The 
MSC/NASTRAN code for this strain energy run is included in Appendix A. 

The results of the strain energy run are shown in Table 3.2. The 
strain energies are given in percent of the total structural strain energy for 
each mode. The "inner vane” elements were the vane elements closest to 
the inner shell. The "outer vane" elements were the vane elements closest 
to the outer shell. The "side outer shell" elements were the outer shell 
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Figure 3.4. MSC/NASTRAN finite element model of test structure. 





Figure 3.5. MSC/NASTRAN finite element 


model of test structure (front view). 


16 


elements closest to the vanes. The data showed that the outer shell was the 
most likely candidate for a damping treatment. The data also showed that 
the vanes were possible secondary candidates for damping treatments. 
The inner shell was confirmed to be an area of very low strain energy and 
was thus determined to be a poor candidate for a damping treatment. Also, 
for both the vanes and the shells, the elements closest to the bolted 
connections gave the highest strain energy. This was interesting since the 


TABLE 3.2. STRAIN ENERGY DISTRIBUTION IN UNDAMPED 
STRUCTURE DETERMINED BY 
MSC/NASTRAN FINITE ELEMENT ANALYSIS. 
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prorat | 102.59 | 121.60 | 26.303 | 132.11] 212.11] 205.03 | 


| AVG/ : 
|mope | 12.82 | 15.20 | 3.29 | 16.51 | 26.51 | 25.63 | 


















experimental mode shapes showed the connections between the vanes and 
the shells to be an area of little deformation. However, since the area near 
the connections is a transition area between an area of little deformation 


and an area of large deformation, large bending stresses are likely (similar 


Lg 


to a cantilever beam near the rigid wall), causing this area to have a high 
strain energy. 


CG FREQUENCY, TEMPERATURE AND MATERIAL TYPES 

The next step in the design process was the determination of the 
applicable temperature and frequency ranges, and the material types for 
the viscoelastic and constraining layers. 

The temperature for all design work was considered to be 70°F. 
Although some small temperature fluctuations occurred during 
experimentation, this temperature was a good average value for use in 
analytical design work. The frequency range of interest was 0 to 300 Hz. 
This low frequency range was chosen since the absorption of sound in 
seawater rises as the square of frequency as shown by the following 
equation [Ref. 9]: 

nab, ie 10° °dB/ yard (3.2) 
where, 
a = absorption of sound in seawater at 65°F 
f = frequency of sound waves 


Thus, sound created by low frequency vibrations is not easily absorbed in 
seawater and must be prevented through the dissipation of structural 
vibration energy. 

The viscoelastic layer material chosen was 3M ISD-112 since it 
displays a very high loss factor in the frequency range of interest. 
Aluminum was chosen for the constraining layer since it is light, relatively 
stiff, easy to fabricate, and is commonly used in most constrained layer 
damping treatments. 


D. DAMPING TREATMENT THICKNESS 

To complete the initial damping treatment design, only the 
thicknesses were left to be determined. To accomplish this, subsystem 
analysis was performed. To do this the structure was broken down into 
smaller subsections, and the damping design was conducted for these 
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smaller subsections. Numerous methods had been suggested for 
accomplishing this subsystem analysis. The system's dimensions could be 
matched to a beam or plate, and then the "effective thickness" could be 
adjusted to make the natural frequencies match [Ref.4 ]. An alternative 
was to insure that the node line distances for the mode of interest matched 
the wavelength of vibration of the beam or plate approximation [Ref. 10]. 
Finally, simply matching the boundary conditions of the subsystem with 
the boundary conditions of a beam or plate had shown reasonable accuracy 
[Ref.11]. 

The test structure's first mode shape, corresponding to the first 
natural frequency of 29.7 Hz is shown in Figure 3.2. From this first mode 
shape, and from close observation of the animated first mode shape, it was 
clear that for the first mode of vibration the outer shell flexed back and forth 
as shown in Figure 3.6. Thus, each quarter of the outer shell behaved 
roughly like the first mode of a curved beam with simply supported ends. 
The equations for the system loss factor for a simply supported beam were 
therefore used to approximately predict system damping for the first mode 
of vibration with one quarter of the outer shell used as the subsystem. 


First Half Cycle of Motion Second Half Cycle of Motion 


Figure 3.6. Motion of test structure during first (29.7 Hz) mode of vibration. 


WS, 


For a three layer system of the type shown in Figure 3.7, which is 
typical of a constrained layer damping treatment, the flexural rigidity can 
be expressed in closed form by: [Ref. 12] 

















ore Et BH, Bly Eee rE eee Di Cee pr 
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H 
31 
EH (H.. - ——]+elE HH. +E.H.H..] 
Dy 22a Zi 2 25 2a 3 3 31 (3.4) 





EoH, 
EH + 9 + g{h H+ EH,+ EH] 
Hi +H, 
it ee (3.5) 
H +2H +H 
= 1 2 3 
A= nee (3.6) 
G, 
o> eS (Sa) 
EHH K 


where, 
H = section thickness 
E = Young's Modulus 
G = shear modulus 
I = area moment of inertia 
K? = wave number 
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and, 
Subscript 1 refers to the base structure 
Subscript 2 refers to the viscoelastic layer 
Subscript 3 refers to the constraining layer 
No subscript refers to the composite structure 


The wave number is different depending on the type of three layer 
structure. For a simply supported beam it is computed as follows: 


Kis ox (3.8) 
where, 
n = mode number 
L = length of beam 





Figure 3.7. Constrained viscoelastic layer damping treatment. 


The preceding equations were used to determine the loss factor and shear 
modulus for the subsystem with a constrained layer damping treatment. 
The following inputs to these equations were used: 


¢ HH, =0.25 in. 

© E, =30X10° PSI 

° Eg =10X10° PSI 

¢ L=Circumference/4 = (1d)/4 = (2*45 in.)/4 = 35.3 in. 


Sn 
¢ Go,N» = found from material property graph 
° Ho,H3 = varied by designer to find optimum values 


The material property values for 3M ISD-112 were displayed on a 
reduced frequency nomograph as shown in Figure 3.8. To find the material 
properties of the viscoelastic material, simply move across the graph from 
the right at the proper frequency, and move down along the constant 
temperature line until the two lines intersect. Then draw a vertical line. 
The intersection of this vertical line with the two material property curves 
gave the desired values for the given frequency and temperature. Thus, the 
values for H2 and H3 were the only variable inputs to the equations for 
system damping and system loss factor. These thickness values were 
varied and loss factor curves for each combination were plotted. A 
FORTRAN program shown in Appendix B was written to solve the system 
loss factor equations for each combination of H2 and H3. The damping 
material thickness, H2, was held constant at 0.005 in. while the 
constraining layer thickness, H3, was varied through eight different values 
from 0.005 in. to 0.040 in., resulting in eight system loss factor curves on 
one plot for a constant H2 of 0.005 in.[Ref. 13]. Then, the damping material 
thickness H2 was changed from 0.005 in. to 0.040 in. in increments of 0.005 
in., with a new graph plotted at each increment. The curves for H2 = 0.005 
in., H2 = 0.030 in. and H2 = 0.040 in. are shown in Figures 3.9 through 3.11. 
Generally, as H3 increased, the loss factor increased for constant H92. Also, 
as H2 increased for a constant H3, the loss factor increased. The loss factor 
curves show that as H2 approached 0.040 in., the curves became very peaky 
and the curves maximum points shifted radically to lower temperatures. 
Also, low values of H2 and H3 did not provide adequate levels of damping. 
Therefore, values of 0.030 in for H2 and 0.040 in for H3 were chosen as final 
design values. At a temperature of 70°F, this combination provided a 
theoretical system loss factor of 0.087. 
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DAMPING PROPERTIES 
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Figure 3.8. Reduced frequency nomograph for 3M ISD-112. 
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Figure 3.9. System loss factor vs. temperature for H9 


frequency is constant at 29.7 Hz. 
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Figure 3.10. System loss factor vs. temperature for f17= 0 Oe0umn, 
frequency is constant at 29.7 Hz. 
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Figure 3.11. System loss factor vs. temperature for Ho = 0.040 in, 
frequency is constant at 29.7 Hz. 
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E. MODAL STRAIN ENERGY ANALYSIS 

To construct the modal strain energy finite element model, the model 
used for the undamped case was modified to include the addition of a 0.030 
in. viscoelastic layer and a 0.040 in. constraining layer to the outer shell. 
The outer shell elements were offset toward the center of the structure by 
0.125 in., and a second layer of grid points was placed 0.030 in. outside the 
first layer. Then, the viscoelastic layer was modeled using CHEXA 
elements, while the constraining layer was modeled using CQUAD4 
elements and offset 0.020 in. away from the center of the structure, as 
shown in Figure 3.12. 


CONSTRAINING LAYER OFFSET 0.020 IN. 


VISCOELASTIC LAYER 


THICKNESS = 0.030 IN. 





BASE STRUCTURE OFFSET - 0.125 IN. 


Figure 3.12. MSC/NASTRAN finite element model 
of constrained layer damping treatment. 
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A modal summation was performed as in [Ref. 2] using the following 
solution to equation (2.7): 


L/w 


OF iP R.i 


X(@) = (3.9) 





where: 


X(@) = amplitude of response as a function of frequency 
QO; = ith natural circular frequency 


Op ; = mode shape value for excitation grid point at ith 


natural frequency 
dp ; = mode shape value for response grid point at jth 


natural frequency 
Nyi= viscoelastic material loss factor 


V..: = strain energy in viscoelastic material at jth 


V,1 
natural frequency 

V;, ; = total strain energy in structure at ith natural 
frequency 

L = magnitude of load applied at excitation grid point 


The loss factor in equation (3.9) is an iterated loss factor. Since it 
changes with frequency over the frequency range of interest, it must be 
adjusted for each natural frequency in the modal summation [Refs. 1,2]. 
The adjustment procedure was carried out as follows. Three finite element 
modal analysis runs were made with three different values of shear 
modulus for the viscoelastic material. This yielded three different natural 
frequencies for each mode, since increasing the shear modulus of the 
viscoelastic material increased the stiffness of the structure. These values 


of shear modulus were added to a plot of shear modulus vs. frequency as 
shown in Figure 3.13. The three shear modulus values were then 
connected, and where the two lines intersected was where the shear 
modulus value was taken. Then, that shear modulus value was used to find 
the new loss factor value from the reduced frequency nomograph, Figure 
3.8. Finally, this loss factor value was adjusted as follows: 


3.10 





where, 


Go REF = shear modulus of viscoelastic material at 


150 Hz 
Go ; = shear modulus of viscoelastic material at jth 


natural frequency 
Nie = loss factor of viscoelastic material from reduced 


frequency nomograph 


This most accurately described the true value of loss factor [Ref. 1]. 
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Figure 3.13. Shear modulus graphical iteration procedure 
for modes 4, 14 and 29. Value for shear modulus is 
chosen at intersection of two curves for each mode. 
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IV. EXPERIMENTS 


A. TESTING ARRANGEMENT 

The testing arrangement is shown in Figure 4.1. The user interface 
was provided by the Hewlett Packard(HP) 3562A Dynamic Signal Analyzer. 
The HP-3562A provided the source signals for the vibration generator, and 
analyzed the incoming input force signals and output acceleration signals. 
The HP-3562A was used to compute the frequency response, coherence, 
input power spectrum and output power spectrum using the discrete 
Fourier Transform at individual frequencies in the swept sine mode. Ten 
averages were used at each data point to provide good coherence 
measurements. Once data taking was complete for a given run, the data 
was saved on disk for later analysis. 

Source signals from the HP-3562A were fed through a Wilcoxon 
Power Amplifier and split by a Wilcoxon Matching Network to power both 
the low frequency magnetic unit and the high frequency piezoelectric unit 
in the Wilcoxon F4/F7 vibration generator. The vibration generator was 
mounted 12 inches from the front and four inches from the inner shell of 
the test structure. The vibration generator was mounted inside a 
waterproof container since earlier research had required the immersion of 
the test structure in water. An integral force transducer was included in 
the base of the vibration generator to measure the input force. This input 
force signal was then fed through a PCB Corporation charge amplifier back 
to the HP-3562A. Two Endevco accelerometers were mounted on the test 
structure to measure system output acceleration. One accelerometer was 
located on the inside of the outer shell of the test structure, four inches from 
the front and twelve inches up from the vane. The other accelerometer was 
located four inches from the front and eight inches from the inner shell on 
the vane of the test structure. One of these output acceleration signals at a 
time could be fed through an Endevco signal conditioner back to the HP- 
3562A. Thus, the HP-3562A provided the source signal, measured input 
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force and output acceleration, computed and displayed the system 
frequency response, and then saved the data for future reference. 


B. TESTING PROCEDURE 
1. Case With Viscoelastic Material on Outer Shell 

The first series of tests was conducted with viscoelastic 
material on the outer shell of the test structure, with one constrained layer 
damping treatment located on each quadrant of the outer shell. Prior to 
application of the constrained layer damping treatment, baseline 
measurements with no damping treatment applied were taken for 
comparison purposes. Then the constrained layer damping treatment was 
applied carefully to prevent the entrapment of air between the viscoelastic 
material and the test structure. The constrained layer damping treatment 
was 18 inches wide, covering the entire width of the shell, and 27 inches 
long, covering most of the area between the area of the bolts as shown in 
Figure 4.2. Seven distinct frequency bands were used, including three wide 
range frequency bands, and four narrow range frequency bands. The three 
wide range frequency bands covered 15 to 115 Hz, 100 to 200 Hz and 200 to 300 
Hz. The four narrow frequency bands covered 25 to 45 Hz, 95 to 115 Hz, 120 to 
140 Hz, and 230 to 270 Hz. The narrow frequency bands were used to "zoom" 
in on natural frequencies of interest to provide higher resolution data at 
these critical areas. Measurements were taken using both the outer shell 
and vane accelerometers. 

2. Case with Viscoelastic Material on Outer Shell and 

Vanes 

The second series of tests was conducted with viscoelastic 
material applied to both the outer shell and the vanes of the test structure. 
Care was again taken when applying the viscoelastic material. The 
damping treatment was again 18 inches wide, covering the entire width of 
the vane, and nine inches long, covering almost the entire distance from 
the inner shell to the outer shell. The same frequency bands were used to 
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facilitate comparison between the baseline, outer shell and outer shell/vane 
cases. 
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Figure 4.1. Testing arrangement. 
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Figure 4.2. Location of viscoelastic material on outer shell. 
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V. RESULTS 


A. FREQUENCY RESPONSE MEASUREMENTS 
1. Case with Viscoelastic Material on Outer Shell 

The frequency responses for the case with viscoelastic material 
on the outer shell are shown in Figures 5.1 through 5.9. In each case the 
dashed curve represents the baseline condition, with the solid line 
representing the test structure after damping treatment application. Also, 
each figure displays amplitude in dB, with 12.5 dB per division on the 
vertical axis, and measurements for Figures 5.1 through 5.7 were taken on 
the outer shell of the test structure. The frequency response of the test 
structure for the 15 to 115 Hz range for both the baseline and viscoelastic 
cases is shown in Figure 5.1. Response amplitudes decreased for all modes 
in this frequency range, with the exception of the modes at 37 and 64 Hz. 
Response amplitude decreases were significant for the modes at 85 and 100 
Hz, with the second mode in each case almost disappearing. Natural 
frequencies shifted lower for all modes with the exception of the modes at 
48, 55 and 64 Hz. The frequency response of the test structure for the 100 to 
200 Hz range for both the baseline and viscoelastic cases is shown in Figure 
5.2. Significant decreases in frequency response amplitude were evident for 
all modes with the exception of the mode at 158 Hz. The mode at 147 Hz 
disappeared. Natural frequencies shifted lower for all modes with the 
exception of the mode at 158 Hz. The frequency response of the test structure 
for the 200 to 300 Hz range for both the baseline and viscoelastic cases is 
shown in Figure 5.3. All frequency response amplitudes for this frequency 
range were significantly lower than the baseline amplitudes. Four modes 
between 240 and 258 Hz disappeared completely. The entire frequency 
response curve shifted down and to the left in this frequency range. A mode 
appeared at 296 Hz that was previously at greater than 300 Hz and thus did 
not appear on the baseline plot. 
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Four narrow range frequency responses were obtained in 
order to increase the resolution at which response amplitude and percent 
damping measurements could be made. The frequency response of the test 
structure for the 25 to 45 Hz range for both the baseline and viscoelastic 
cases is shown in Figure 5.4. In this range, two out of three of the modes 
decreased in response amplitude, and all modes shifted lower in natural 
frequency. The natural frequency for the mode at 29.6 Hz dropped to 29.0 
Hz, and its response amplitude dropped 4.6 dB. The natural frequency for 
the mode at 35.8 hz dropped to 35.6 Hz, and its response amplitude 
increased by 7.9 dB. The frequency response of the test structure for the 95 to 
115 Hz range for both the baseline and viscoelastic cases is shown in Figure 
5.5. Both modes in this range showed significant natural frequency shifts 
and amplitude response decreases. The natural frequency of the mode at 
99.4 Hz shifted to 97.4 Hz, with the response amplitude decreasing by 16.0 
dB. The natural frequency of the mode at 100.8 Hz shifted to 98.85 Hz, with 
the response amplitude decreasing by 12.6 dB. The frequency response of 
the test structure for the 120 to 140 Hz range for both the baseline and 
viscoelastic cases is shown in Figure 5.6. The double mode of the baseline 
structure shows up as one mode in the damped structure. The natural 
frequency for the mode at 126.3 Hz shifted to 123.8 Hz, with the response 
amplitude dropping by 12.4 dB. The natural frequency of the mode at 136.2 
Hz shifted to 133.1 Hz, with the response amplitude dropping by 17.8 dB. 
The frequency response of the test structure for the 230 to 270 Hz range for 
both the baseline and viscoelastic cases is shown in Figure 5.7. The 
response amplitude of every mode in this range has been significantly 
reduced, with four of the six modes in this range disappearing after the 
application of the damping treatment. The natural frequency of the mode at 
264 Hz decreased to 261.2 Hz, with the response amplitude decreasing by 
14.5 dB. 

The results for the measurements taken on the vane of the 
structure are generally the same as for those taken on the outer shell of the 
structure. As an example, the frequency response for the 15 to 115 Hz run 
is shown in Figure 5.8. As was the case in Figure 5.1, the higher frequency 
modes at 84 and 100 Hz showed large drops in response amplitude, while 
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the drops in response amplitude for the low frequency modes at 29 and 36 
Hz were small. For the narrow frequency range run from 230 to 270 Hz 
shown in Figure 5.9, the vane measurement was similar to the outer shell 
measurement shown in Figure 5.7. All modes either disappeared or had 
their response amplitudes greatly reduced. Even the mode at 243 Hz, which 
shows up strongly in the vane measurement, was reduced in amplitude by 
19 dB. 
2. Case with Viscoelastic Material on Outer Shell and 

Vanes 

The frequency responses for the case with viscoelastic material 
on the outer shell and on the vane are shown in Figures 5.10 through 5.18. 
The curves were different from those taken with viscoelastic material on 
the outer shell only since modal response amplitude reductions and 
increases were magnified slightly by the addition of the viscoelastic 
material on the vanes. For instance, the response amplitude of the mode at 
100 Hz, as seen in Figure 5.14, decreased by 17.5 dB versus 16 dB for the case 
with viscoelastic material only on the outer shell. Also, the response 
amplitude for the mode at 264 Hz, as shown in Figure 5.16, decreased by 
16.5 dB versus 14.5 dB for the case with viscoelastic material only on the 
outer shell. Interestingly, the 36 Hz mode that experienced a sizable 
increase in response amplitude due to the application of the damping 
treatment on the outer shell had its response amplitude increased even 
more for the case of viscoelastic material on both outer shell and vane. The 
response amplitude for this case, as seen in Figure 5.13, increased 12.5 dB 
versus 7.9 dB for the case with viscoelastic material on the outer shell only. 

Measurements taken on the vane for this case were similar to 
the measurements taken on the outer shell. Response amplitudes generally 
decreased similar amounts and frequency shifts were also similar. 


B. COMPARISON WITH MODAL STRAIN ENERGY METHOD 
The frequency response results for the Modal Strain Energy(MSE) 

Method versus the case of viscoelastic material on the outer shell only are 

found in Figures 5.19 through 5.25. The correlation between the frequency 


response predicted by the MSE method and the results obtained 
experimentally was not good. Natural frequencies predicted by the MSE 
method were lower than those obtained experimentally for frequencies less 
than the mid range (150 Hz) frequency, and natural frequencies predicted 
by the MSE method were higher than those obtained experimentally for 
higher than mid range frequencies. The MSE curves generally showed 
amplitudes higher than those found by experimentation, and generally 
showed less damping. Also, occasional modes appear in the MSE curves 
that did not show up in testing. Note that the experimental curve in Figure 
5.23 actually has its mode at 97 Hz, and is thus approximately 10 Hz higher 
than the MSE prediction. 


eo STRUCTURAL DAMPING 

Structural damping was measured for four modes for each damping 
treatment case using the half power bandwidth method. Structural 
damping as a percentage of critical damping is shown in Figure 5.26. All 
cases are for measurements taken on the outer shell. As shown in Figure 
5.26, the undamped test structure is an extremely lightly damped 
structure. Addition of the constrained layer viscoelastic damping treatment 
to the outer shell increased percent damping by 1.5 to 8 times. Further 
addition of the damping material to the vane resulted in marginally higher 
damping in three out of four cases.The highest damping recorded was 1.044 
percent for the mode at 123.7 Hz with damping material applied to the outer 
shell and vane. Many modes, especially in the 230 to 270 Hz range, could not 
be plotted since addition of damping treatments made the modes completely 
disappear or made them rise less than the 3 dB required to get a damping 
measurement using the half power bandwidth method. The MSE method 
underpredicted the amount of damping obtained. 


D. DISCUSSIONS 

The addition of a constrained layer viscoelastic damping treatment 
was seen to be a very effective method of controlling large resonant 
vibrations in the test structure. Most modes had response amplitude 
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reductions in the 15 dB area, with some modes disappearing completely. 
This was accomplished with a damping treatment which weighed 12.9 
pounds, or approximately two percent of the total weight of the 
approximately 650 pound structure. Also, structural integrity of the test 
structure was not adversely affected by the addition of the damping 
treatment. 

Addition of the damping treatment to the vane as well as the outer 
shell was found to be only marginally more effective in controlling resonant 
vibrations. Response amplitude reductions were only slightly larger with 
the addition of the damping treatment to the vane. This was expected due to 
the large percent of test structure motion that occurred in the outer shell 
area, as determined by the mode shape and strain energy analyses. In 
summary, if most of the motion occurs in the outer shell, most of the energy 
dissipation will occur there also. 

Vibration energy dissipation was worst in the low frequency region 
below 75 Hz. This was due to the low shear modulus of 3M ISD-112 in this 
area. As equation 2.11 shows, the energy dissipation in a constrained layer 
system is directly dependent on the shear modulus. As Figure 3.13 shows, 
the shear modulus of 3M ISD-112 in the area from 20 to 75 Hz is only 50 to 
100 PSI. This is simply too low to provide good energy dissipation. In effect, 
the viscoelastic material just "gives", causing great reductions in possible 
energy dissipation. 

The response amplitude for the second mode at 36 Hz actually 
increased with the addition of the constrained layer damping treatment. As 
seen in Figure 3.3, the outer shell did move in this mode but experienced 
almost no deformation, while the vane experienced a large bending 
deformation. Thus, application of a damping treatment to the outer shell 
provided no energy dissipation at all, and the added mass of the damping 
treatment caused the vane to experience a larger deformation. This caused 
the response amplitude to increase. Likewise, addition of the damping 
treatment to the vane only increased the mass more, causing any energy 
dissipation which may have occurred in the vane to be overshadowed by the 
increases in response amplitude. Also possibly contributing was the 
presence of entrapped air in both the outer shell and vane damping 
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treatments. These air pockets contributed to a reduction in damping 
treatment effectiveness [Refs. 2,14]. Although difficult to quantify, the 
presence of imperfections in the viscoelastic bond to the test structure could 
have contributed to the increases in response amplitude seen in the second 
mode. Of course, energy dissipation at all modes would probably improve if 
the air pockets could be completely removed. 

The MSE method did not give a good prediction of the experimental 
results obtained. In general, the MSE method tended to show response 
amplitudes that were too high with too little damping. In addition, some 
extra modes appeared in the MSE method results that did not appear 
during experimentation. The reasons for these results are as follows: 


¢ The bolts which connect the inner and outer shells to the vanes of the 
test structure were not modeled in the finite element analysis. The 
finite element analysis assumed that the test structure was a 
smooth, continuous structure. Actually, the structure contained 144 
bolts, with 96 bolts contained in the inner shell, and 48 bolts contained 
in the outer shell. These bolts acted as stiff clamps, holding the vanes 
and outer shell in place, thereby reducing the response amplitudes of 
those sections of the test structure. Since the finite element model did 
not model these bolts, the MSE method predicted higher responses 
than were obtained. 

e The aspect ratio is the ratio of the length of the longest to the shortest 
side for a finite element. An aspect ratio of one is the most desirable 
case since high aspect ratio elements cause large errors to be 
introduced into a finite element analysis. The aspect ratio for the 
CHEXA elements used to model the viscoelastic material of the test 
structure was 294.7. Although Kienholz and Johnson have reported 
adequate results for simple structures with aspect ratios this high 
[Ref. 1], the high aspect ratio warning message issued by 
MSC/NASTRAN during the analysis indicates that this was a source 
of error. 

¢ When modeling a curved structure using flat elements, the included 
angle of the flat elements should be 10 degrees or less to obtain 
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displacement errors of four percent or less [Ref. 15]. In fact, even 
smaller included angles may be required when performing modal 
analysis as was performed for this analysis [Ref. 15]. The included 
angle for the outer shell elements of the finite element model was 22.5 
degrees, and the included angle for the inner shell elements was 45 
degrees. This coarseness in the finite element model undoubtedly 
contributed to the less than exact results obtained with the MSE 
method. 

° The viscoelastic material properties used in the finite element model 
were found from Figure 3.8, which was provided by 3M Corporation. 
However, the material property curves used by Maurer [Ref. 2], found 
in Vibration Damping [Ref. 4] and used by Drake [Ref.13] in 
designing damping treatments with 3M ISD - 112 are all different. In 
addition, Drake and Kluesener [Ref. 16] showed that different 
material property curves were obtained depending on how the 
properties were measured. In summary, there is no consensus on 
exactly what the material properties of 3M ISD-112 are, and thus it is 
possible that the material property values used in the finite element 
model were inaccurate. This possibility is reinforced by the fact that 
the shear modulus found from Figure 3.8 was low compared to the 
value found from the other sources stated above, and the damping 
predicted by the MSE method turned out to be low also. 

e Viscoelastic material property values of different samples of the 
same material tend to vary slightly. Thus, even if the finite element 
model used material property values that were accurate as average 
values, it is possible that the actual viscoelastic material used during 
experimentation had material property values that differed from the 
values used in the prediction. 


The results obtained by the MSE method were conservative. Response 


amplitudes were generally higher than those obtained by testing, and 
damping was less than testing showed. 
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Figure 5.1. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell area, 15 to 115 Hz. 
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Figure 5.2. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell area, 100 to 200 Hz. 
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Figure 5.3. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell area, 200 to 300 Hz. 
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Figure 5.4. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell area, 25 to 45 Hz. 
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Figure 5.5. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell area, 95 to 115 Hz. 
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Figure 5.6. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell area, 120 to 140 Hz. 
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Figure 5.7. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell area, 230 to 270 Hz. 
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Figure 5.8. Frequency response measured on vane of test structure, 
damping treatment applied to outer shell area, 15 to115 Hz. 
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Figure 5.9. Frequency response measured on vane of test structure, 
damping treatment applied to outer shell area, 230 to 270 Hz. 
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Figure 5.10. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell and vane areas, 15 to115 Hz. 


o2 






ea 


Pp] tty 
| be? SERRE 





Damping Treatment on Outer Shell and Vane 


eerneeceeeee---- [Jndamped (Baseline) Test Structure 


N 
Ble 
Oo 
U) 
Lu 
Y ) 
(a) 
(3 - 
Lu 
ie @ Wl SI 
ie, . oA LN) 
7) le) a U 
= = ON Mm | xX 
aie U lif 


Figure 5.11. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell and vane areas, 100 to 200 Hz. 
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Figure 5.12. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell and vane areas, 200 to 300 Hz. 
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Figure 5.13. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell and vane areas, 25 to 45 Hz. 
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Figure 5.14. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell and vane areas, 95 to 115 Hz. 
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Figure 5.15. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell and vane areas, 120 to 140 Hz. 
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Figure 5.16. Frequency response measured on outer shell of test structure, 
damping treatment applied to outer shell and vane areas, 230 to 270 Hz. 
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Figure 5.17. Frequency response measured on vane of test structure, 
damping treatment applied to outer shell and vane areas, 15 to115 Hz. 
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Figure 5.18. Frequency response measured on vane of test structure, 
damping treatment applied to outer shell and vane areas, 230 to 270 Hz. 
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Figure 5.19. Modal Strain Energy method vs experimental frequency 
response, damping treatment applied to outer shell area, 15 to 115 Hz. 
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Figure 5.20. Modal Strain Energy method vs experimental frequency 
response, damping treatment applied to outer shell area, 100 to 200 Hz. 
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Figure 5.21. Modal Strain Energy method vs experimental frequency 
response, damping treatment applied to outer shell area, 200 to 300 Hz. 
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Figure 5.22. Modal Strain Energy method vs experimental frequency 
response, damping treatment applied to outer shell area, 25 to 45 Hz. 
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Figure 5.23. Modal Strain Energy method vs experimental frequency 
response, damping treatment applied to outer shell area, 95 to 115 Hz. 
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Figure 5.24. Modal Strain Energy method vs experimental frequency 
response, damping treatment applied to outer shell area, 120 to 140 Hz. 
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Figure 5.25. Modal Strain Energy method vs experimental frequency 
response, damping treatment applied to outer shell area, 230 to 270 Hz. 
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Figure 5.26. Structural damp 


VI. CONCLUSIONS 


The addition of a constrained viscoelastic layer damping treatment to 
a large, lightly damped structure was found to be very effective in reducing 
unwanted structural vibrations. The large reduction of response 
magnitude was encouraging, especially considering the light weight of the 
applied damping treatment. The relatively poor performance of the 
damping treatment in the lowest fifth of the frequency range studied 
indicated that viscoelastic materials with low shear moduli in certain 
frequency ranges are not effective in absorbing vibrational energy. 

The use of visual mode shape observation and strain energy analysis 
to determine suitable locations for damping treatments was found to be 
effective. These analyses predicted that the most damping would be 
obtained by placing a damping treatment on the outer shell and 
experimentation proved that prediction to be true. In addition, when the use 
of these analyses indicates that an area of a structure is a poor candidate 
for a damping treatment, simply adding more damping material to that 
area may provide only marginally more damping, and may not be effective 
in light of the added mass, expense, and time required for application. 

The Modal Strain Energy method is a viable means of predicting the 
frequency response of a structure with an applied damping treatment. 
However, the finite element model itself must be refined enough and must 
be accurate enough to allow the MSE method to perform adequately. 
Indeed, the need to accurately model complex structures makes the MSE 
method all the more attractive since the MSE method requires only three 
modal analysis runs, while direct frequency response techniques may 
require many more runs than that to adequately describe the behavior of the 
structure. In addition, since frequency response analysis takes much more 
time than modal analysis, direct frequency response techniques would not 
be cost effective for many complex structures. 
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VIT. RECOMMENDATIONS 


The test structure used in this investigation was very much like 
many modern structures in that it was large, complex and lightly damped. 
Thus, it is recommended that this structure be used as a test unit for 
further damping research. For instance, it has been suggested that 
concentrated dampers, similar to but smaller than the ones used on large 
buildings, be installed between the inner and outer shells of the test 
structure to take advantage of the relative motion between these two parts of 
the structure.* Many other innovative damping treatments are possible for 
use on this test structure, and the instrumentation is in place to quickly 
record their effects. 

A more refined and accurate finite element model of the test 
structure should be developed that includes the effects of the bolts. This 
could include the use of direct matrix input to model the added stiffness and 
mass of the bolts. Also, the use of singly curved shell elements could be 
used to more accurately model the inner and outer shells. Finally, an 
alternate method of modelling the viscoelastic material such as using beam 
elements instead of solid elements should be attempted. The MSE method 
can then be retested using that refined model for comparison to the 
experimental data obtained. 

The use of multi-layered constrained layer damping treatments 
should be investigated. Multi-layered damping treatments can increase the 


* Conversation with Dr. Dwayne Nelson, 3M Corp., Minneapolis- 
St.Paul, MN, 8 Aug 1988. 
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amount of damping that can be obtained, and they can increase the 
conformability of the treatment since the constraining layers can be made 
of a less stiff material [Ref. 13]. The increased conformability of the 
treatment may make application easier and may therefore result in fewer 
bond imperfections between the treatment and the base structure. 
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APPENDIX A 


INPUT DATA SETS FOR MSC/NASTRAN FOR 
UNDAMPED STRAIN ENERGY ANALYSIS AND FOR 
DAMPED MODAL STRAIN ENERGY(MSE) ANALYSIS 


The major difference between the two sets of NASTRAN input data is 
that the second set includes another set of grid points above the outer shell. 
In between this second set of grid points and the original grid points are the 
solid CHEXA elements used to model the viscoelastic material. Also, in the 
second case the CQUAD4 elements in the outer shell and constraining 
layer are offset from the grid points to make room for the viscoelastic 
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Figure A.1. MSC/NASTRAN data set for 
undamped strain energy run (Page 1). 


73 


ly 


shaoe 


AP ids 55, lr te?) ones 
=, (V jp 237 CS 
=(6) 

Arid, bsp Volo onl, renee 
=,*(1),=,=,*(99.),= 
=(?) 

Qridvy 7 Up lw U 3s ete, Cece 
=,*(1),=,=,*(90.),= 
=(?) 

aridg O@Vel, ce soos lee 
=,*(1),=,=,*02°2.5),5 

= (oe 

fo hae Wik ae ie ae ee 
Sk Cl )ip =p ae lee 
=(6) 

qgria, 105,11, 1 1a) 0450 eee 


=,e(1),57=-4(99.),5 

=(2) 

Jerid, 19901, 12 le4%e7 Cee 
=,*(1),=7=9 *19%% 7 = 

=(?) 

Seidel |! See © ou ee emma 
Sh ea = ee 

=(14) 

equad4,1-1971,335, 34-7 

=, OC) )e tet (let Cli tly ew 
(>) 

hel Op ae GOs 5359! 
Prarie A Se Asean slern. S 

seC1 )p Sr © Cle fle st 1 
Cso 

pe Nie 10a Ui Coto oer 

pT NOS ee oe oO 

= el) pap et ype ae ee Ls el 
=(5) 

Sere Oe Chet oS 
rey ee On tert ree: 

=k (1 jopay tl 2) (1 ge Cee 
=(?) 

STP open el oe tj) 

= pe (Vy 80 1 ee On ae 
=(2) 

S755 Oe Poel ie teu 

= e(1),=,2 01) * (204 sew 
=(?) 


Figure A.2. MSC/NASTRAN data set for 
undamped strain energy run (Page 2). 
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Figure A.3. MSC/NASTRAN data set for 
undamped strain energy run (Page 3). 
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Figure A.4. MSC/NASTRAN data set for damped 
Modal Strain Energy analysis (Page 1). 
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Figure A.5. MSC/NASTRAN data set for damped 
Modal Strain Energy analysis (Page 2). 
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Figure A.6. MSC/NASTRAN data set for damped 
Modal Strain Energy analysis (Page 3). 
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Figure A.7. MSC/NASTRAN data set for damped 
Modal Strain Energy analysis (Page 4). 
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Figure A.8. MSC/NASTRAN data set for damped 
Modal Strain Energy analysis (Page 5). 
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APPENDIX B 


FORTRAN PROGRAM USED TO COMPUTE LOSS 
FACTORS FOR THREE LAYER SYSTEM 


Calculations of loss factors for the three layer subsystem used to 
determine damping material thicknesses were performed using a 
FORTRAN prorgam. The program used to perform these calculations is 
shown in this appendix. The program reads in values from a data table of 
shear modulus and loss factor versus temperature for 3M ISD-112, and 
outputs the loss factor of the composite system. Each run of the program 
supplied data of loss factor versus temperature for eight values of 
constraining layer thickness. These values would then be plotted as shown 
in Figures 3.9 through 3.11, for different values of viscoelastic material 
thickness. 
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PROGRAM CARPET 
DE IE IE HEHEHE HEHE HE HE HE IE DE HE 3E 3 HE DE IED 3 3 HE HE HE DE DE DE DE 3 3 IE IE IE DE DE DE IE BEBE DE IE DE 3 IE 3E 3 IE IE DK IE DE FE IE HE HC HE HE HE HE DE HE EE DEE IE 
® PROGRAM CARPET PRODUCED BY LT. JAMES R. NAULT USN . THIS PROGRAM xX 
% READS IN THE VALUES OF LOSS FACTOR AND SHEAR MODULUS FOR ISD-112 x 
® AND USES THESE VALUES TO COMPUTE A CARPET PLOT FOR THE SYSTEM OF~ x 


% INTEREST. x 
HE HE 3 3 3 HE HE 3 HE BE BE HE DE DE BE HE HE DE DE IE HE 3 DE HE HE IE HE HE DE DE DE DE DE DE DE DE DE DE DE HE HE IE FE IE IEE DE 3E 3 IE DE DE HE IE DE DE DK HE HE HE IE HE IE HE DE 3 IE HE HE 


REAL F,FROM,MROM,N,ML,ETAFRO,SL,SH,FROL,C,T0,TEMP,A,LOGFR,FR 
REAL LOGM,M,LOGETA, ETA2,LOGET1,H1,E1,H2,H3,E3,K,G,LITA,LITB 
REAL ETAS(8),EHCUBE,E,H,H31,AREL,AIMG,PI 
® ASSIGN FREQUENCY OF INTEREST VALUE 
F=2977 0 
* COMPUTE SEs 
PI=4.0xXATAN(1.0) 
® ASSIGN CONSTANT VALUES FOR STRUCTURE 
H1=0.25 
E1=30000000.0 
E3=10000000.0 
K=Fi7S5-3 
H2=0.030 
* READ IN MODULUS AND LOSS FACTOR FOR ISD-112 
DO 10 J=0,100 
READX, TEMP,M,ETA2 
*% COMPUTE SYSTEM YOUNG'S MODULUS AND SYSTEM LOSS FACTOR 
DO 20 I[=1,8 
H3=REALCI)7200.0 
H=H1+H2+H3 
H21=CH1L+H2)72.0 
Hol=Chl+2 OAnicti 7 ac 
G=M/7 CES XH3S XHZEKEX2 ) 
E2=3.0XM 
LITA=E1*H1+G¥CE1¥H1+E3XH3) 
LITB=ETA2XGXCE1XHI+E3X%H3) 
AREL=GXE1 XH1¥XE3XH3XH31¥X2HCLITA+ETACHLITB) 
AIMG=GXE1*H1 XESXH3¥H31¥*2*C ETAZXLITAD 
BREL=E1XHIXE2XH2XHSLHCLITA+LITBXETA2) 
BIMG=E1 XH1¥XE2RH2XHS1XCLITAXETA2-LITB) 
CREL =2XGREZRH2XE3XH3SHH21] XHSLHCLITAXC1l-ETAZ#X2)+24LITB 


+ XETAZ)- 
CIMG=2*XGREZXH2XESXH3XH2Z1XH3S1 XC 2¥LITAXETAZ2-LITB+LITB® 
+ ETA2**2) 
EHCUBE=E1] XH] ¥¥3+E3¥H3¥¥3412.0%(AREL-BREL-CRELI/ CLITA 
+ X¥2+LITB¥¥2) 


E=EHCUBE/ (HXX3) 
ETAS(CIJ=12.0¥CAIMG-BIMG-CIMG)/CCLITAXX¥2+LITB*¥2) 
+ XEHCUBE) 
20 CONTINUE 
PRINT 15, TEMP, ETASC1), ETAS(Z), EIAS CS EIAoc a 
+, ETASCS)EIASCO)FEIASCZ LIAS GoD 
5 FORMATC1X,F6.1,2X,8F09.5) 
0 CONTINUE 


END 


pms ps 


Figure B.1. FORTRAN program to calculate system loss factor. 
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